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Abstract 

We prove that the Faddeev-Popov ghost dressing function in the Yang- 
Mills theory is non-zero and finite in the limit of vanishing momenta and hence 
the ghost propagator behaves like free in the deep infrared regime, within the 
Gribov-Zwanziger framework of the D-dimensional SU(N) Yang-Mills theory 
in the Landau gauge for any D > 2. This result implies that the Kugo-Ojima 
color confinement criterion is not satisfied in its original form. We point out 
that the result crucially depends on the explicit form of the non-local horizon 
term adopted. The original Gribov prediction in the Landau gauge should be 
reconsidered in connection with color confinement. 
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1 Introduction and main results 



The behaviors of the gluon and ghost propagators in the deep infrared regime is im- 
portant for clarifying the dynamics of non-perturbative nature in non-Abelian gauge 
theories. In particular, it is believed that they are intimately related to quark, gluon 
and more generally color confinement [HE]- 

From the viewpoint of the analytical treatment of the quantum field theory, the 
exact path integration of the quantum Yang-Mills theory in the continuum spacetime 
is still difficult in the sense that the fundamental modular region is not yet specified 
in a manageable form of the Lagrangian or Hamiltonian field theory in the Lorentz 
covariant gauge. Nevertheless, it will be important to consider how to incorporate the 
existence of the Gribov horizon as the boundary of the 1st Gribov region to remove 
the Gribov copies as much as possible 0,IU[5j[6] • 

In a previous paper [7] , we have studied an effect on the ghost propagator or ghost 
dressing function G(k 2 ) of the horizon condition which plays the role of restricting 
the functional integral region of the gluon field to the 1st Gribov region. We have 
rewritten the Zwanziger horizon condition in terms of the ghost dressing function and 
the Kugo-Ojima parameter u(0) for color confinement. This has enabled one to study 
which value of the Kugo-Ojima parameter u(0) is allowed if the horizon condition 
is imposed. Although all the calculations were performed in the limit of vanishing 
Gribov parameter 7 for simplicity, the obtained value is consistent with the result 
of numerical simulations. In fact, the direct measurements on a lattice [HIE] show 
u(0) = —0.6 ~ —0.8. Consequently, the ghost propagator behaves like free and the 
gluon propagator is non-vanishing at low momenta. 

In this paper, we study the same issue in the D-dimensional SU(N) Yang-Mills 
theory in the Landau gauge within the Gribov-Zwanziger framework with a non-zero 
Gribov parameter 7 7^ 0. We discuss how the restriction of the integration region to 
the (1st) Gribov region constrains the possible value for the Kugo-Ojima parameter 
for color confinement and the ghost dressing function. 

We prove that the ghost dressing function G(k 2 ) is non-zero and finite in the 
deep IR limit k — > irrespective of the number of color N, and hence the ghost 
propagator behaves like free in the deep infrared regime for D > 2. In addition, with 
an additional input, we have 

G{k 2 ^ 0)^3 (D = 4). (1.1) 

This result is equivalent to say that the Kugo-Ojima color confinement criterion 
u(0) = — 1 is not satisfied in its naive form. Rather, we find the exact value for 
the KO parameter 

u(0) = -2/3 (D = 4), (1.2) 

irrespective of the number of color. These results are in harmony with the decoupling 
solution of the Schwinger- Dyson equation [IT] , recent lattice results [T2l[T3l[T^[To'l[T6"l 
[T7] . and other approaches [T8lfT9l [20] . 
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2 General settings 



We consider the Gribov-Zwanziger theory [I] for the D- dimensional Euclidean SU(N) 
Yang-Mills theory in the Landau gauge, which is defined by the partition function: 

Z 7 := J [ds/\8{& t s^ li ) detKexp{-S YM -J J d D xh(x)}, (2.1) 

where Sym is the Yang-Mills action, 

S Y mW\ ■= J d D x-^ u ■ J^, ^ := dps** - dvsfp + x st v , (2.2) 

K is the Faddeev-Popov operator associated with the Landau gauge fixing condition 
<9^(x) = 0: 

K\sf\ := -d^D^} = -d^ + g^x), (2.3) 
h(x) = h[e/](x) is the Zwanziger horizon function given by0 

h(x) := J d D ygf ABC ^ B (x)(K- 1 ) CE (x,y)gf AFE ^ F (y). (2.4) 

In addition, the parameter 7 called the Gribov parameter is determined by solving a 
gap equation, commonly called the horizon condition: 

(h(x)) 1 = (N 2 - 1)D. (2.5) 

Here the dot and the cross are defined as 

srf ■ 38 := s^ A ^ A , x m) A := f ABC #/ B & c , (2.6) 

using the structure constant of the gauge group G = SU(N). 

It should be remarked that the action corresponding to the partition function (12.1 p 
contains the non-local horizon term: 

J d D xh{x) := J d D x J d D ygf ABC ^ B {x){K~ 1 ) CE {x,y)gf AFE ^{y). (2.7) 

In what follows, we assume that the Euclidean version can be continued analytically 
to the Minkowski version by the Wick rotation. 

Moreover, it has been shown p3,[6] that the non-local action ( 12. ip can be put in an 
equivalent local form by introducing a set of complex conjugate commuting variables 
£, £ and anticommuting ones u, u>, which we call the localized Gribov-Zwanziger (GZ) 
action. In fact, the horizon term can be rewritten into the local form by introducing 
auxiliary fields (See e.g., [2"T]): 

e -7/«Prt(«) = J [dfl^fafidQ] exp {-S,[^,Z,luj,u}} , (2.8) 

1 We have changed the convention. In [7], 7 must be replaced by —7 in (1.1) and (2.1). In |7J, 
the minus sign in (1.2), (1.4), (3.1a), (3.1b) and (3.1c) should be removed. The minus sign should 
be inserted in (3.2). After these changes, the results arc unchanged. 
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where 



S^^luM =: J d D x[^ A K AB $ 



- uj ca K ab lv cb 



+ i^ 9 f ABC ^tf + n Xl2 9f ABC ^tf\ (2-9) 
Thus, the localized Gribov-Zwanziger theory is given by 

Z GZ ■= J \dst\ [dSS\ [cfif] [dtf] [d£] [d£] [duo] [du>] exp{-% - 5 7 }, (2.10) 

where we have introduced the Nakanishi-Lautrup auxiliary field 3$^ the Faddeev- 
Popov ghost field ^ and the antighost field ^ , 

S G z :=^K tf, <*f, 38\ + £ 7 K, e, I u, ou] 

=5 YM K] + S GF+FP [^, V,<e,&\ + S y [j*,Z,Z,u,iD], (2.11) 
with the gauge fixing and the Faddeev-Popov term of the BRST exact form: 

Sgf+fp := / d D x {SS ■ + i<i ■ dpDje) = -S [itf ■ (d^)] , (2.12) 

using the BRST transformation 8. The localized GZ theory is renormalizable to all 
orders of perturbation theory. Hence, the restriction to the (first) Gribov region Q 
makes perfect sense at the quantum level, and finite results are obtained consistent 
with the renormalization group. In the limit 7 — > 0, the integrations over the auxiliary 
fields cancel and the GZ theory reduces to the usual BRST approach for the Yang- 
Mills theory: 

Z YM ■= J [d£/] [d38\ [<fif] [d&\ expi-Sph}. (2. 13) 

It is known that the localized horizon term breaks the usual BRST invariance of 
the GZ theory. Therefore, the GZ theory is no longer invariant under the usual BRST 
transformation. In other words, the BRST symmetry is not the exact symmetry for 
the GZ theory (unless 7 = 0). Nevertheless, one can find the modified "BRST" 
symmetry for the GZ theory with nilpotency [21] (or without nilpotency [22]), which 
reduces to the usual BRST symmetry in the limit 7 — > 0. However, the modified 
"BRST" symmetry inevitably becomes non-local even in the localized GZ theory. 

It should be noted that the Kugo-Ojima (KO) color confinement criterion was 
obtained in the framework of the usual BRST quantization for the Faddeev-Popov 
approach, which corresponds to the 7 = case of the above Gribov-Zwanziger for- 
mulation. The usual BRST approach does not take care of the Gribov copy problem. 
Therefore, if one begins to avoid the Gribov copy by restricting the space of gauge 
field configuration, it may happen that the Kugo-Ojima criterion u(0) = — 1 based 
on the usual BRST approach does not necessarily hold. 



3 Horizon function and ghost propagator 

In what follows, we define the Fourier transform of the two-point function for com- 
posite operators by 

<#Vf > fc := / d°xe^)(O|T[^(x)0f (y)]|0). (3.1) 
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In a previous paper [7j, the 7 = case has been studied and the following identity 
has been derived. The average of the horizon function (assuming the translational 
invariance) is exactly rewritten in terms of three parameters G(0), u(0) and w(0) 
defined below: 



(h(0)) =V D J d u x{h{x)) 

= -(N 2 - 1) {(D - l)u(O) + G(0)[u(0) + w(0)}} , (3.2) 



where we have defined the volume Vo of the D-dimensional Euclidean space (time). 
This expression given in a previous paper is same in the bare case before the renor- 
malization as 

(h(0)) = -(N 2 - 1) {Du(0) + w(0) - G(0)[u(0) + w(0)] 2 } . (3.3) 

Here G(0) is the IR limit k 2 — > of the ghost dressing function G(k 2 ) defined from 
the ghost propagator l^€ Ac € B ) k by 

G(k 2 )5 AB := -k 2 {^ A( i B ) k , (3.4) 

u(0) is the Kugo-Ojima parameter defined by the k 2 — > limit of the Kugo-Ojima 
function u{k 2 ): 

{{D^) A {gK x V) B ) k := L u - k -^f\ 5 AB u(k 2 ), (3.5) 



and u>(0) is the massless pole residue in 

((g^ x V) A {gstf v x tf) B )f PI 



irAB 



(3.6) 



where the modified one-particle irreducible (mlPI) part of ((g^ x ( ^') A (g^' l , x c &) B ) k 
is defined by (the diagrammatic representation is give later.) 

((g^ x ^) A {gK x ^) B ) k 

■■=((9< x V) A (9< x V) B )f pl + ((^ x ^) A ^ c ) k m {^ D ) k {^ D (gK x <£) B )Y\ 

(3.7) 

with the IPI part of {(g^ x ^) Ac i B ) k and (^(g^ x ^) B ) fc defined by 

x ^)V B ) fe x ^tf*)*, (3.8) 

(<^(<^ x ^) B ) fe H^^HO^GM x ^) B )i pi - (3.9) 

In the Landau gauge, especially, the ghost dressing function G{k 2 ) is related to 
the two functions u(k 2 ) and w(k 2 ) as 

G(k 2 ) = [l + u(k 2 )+w(k 2 )}- 1 . (3.10) 

In the Landau gauge, there is a symmetry for the exchange between ghost and 
antighost, called the Faddeev-Popov conjugation invariance. See Appendix A. 
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In the next section, we confirm that all the relations in the above hold also in 
the 7 7^ case of the Gribov-Zwanziger theory, although they are originally derived 
in the 7 = case [7j. Therefore, the conclusions derived from them in the previous 
paper [7] are valid also in the 7 ^ case of the Gribov-Zwanziger theory. Although 
the Gribov-Zwanziger theory involves 7, these relations do not involve the explicitly 
7-dependent extra terms and the 7 dependence appears only through the functions 
F, u and w implicitly. 

From (E2D or lEOl . due to the horizon condition (h(0)) = (N 2 - l)D < +00, it 
is obvious that G(0) is finite as far as u(0) and w(0) are finite. The ghost dressing 
constant obeys 

(D - IMP) + (h(0))/(N 2 - 1) (D — IMP) + D 
1 ' u(0)+w(0) u(0)+w(0) ' 1 ' ' 

On the other hand, we have from (13. 101) 

G(0) = [l + u(0) + w(0)}-\ (3.12) 

If u(0) + w(0) = P, then G(0) is finite, i.e., G(P) = 1, due to f[3~T2D . In this case, 
we have u(0) = — w(0) = —Dj(D — 1) from the consistency with (13. lip . If u(0) + 
w(0) 7^ P, then G(0) can not be divergent and is finite due to the horizon condition 
(h(0)) = (N 2 — 1)D < +00, see (13. 1 II) . In any case, thus, the ghost dressing function 
G(k 2 ) is finite in the deep IR limit. 

It was implicitly assumed in Kugo [2] that H 

iw(0) = 0. (3.13) 

It should be checked whether w(0) = is true or not. According to numerical 
simulations on a lattice [9l[TP] (see Figure 5.4 in section 5.2.2 [9]) and an independent 
study [23] based on [SUES], w(0) = seems to be true. In this case, 

G(P) = -D + 1 + (3.14) 

In the case of w(0) = P, therefore, we have two relations: 

• a relationship (h(0)) = -(N 2 - 1) \-Du(0) + U ^ . 1 , (3.15a) 

{ 1 + "(0) J 

• the horizon condition (h(0)) = {N 2 - 1)D. (3.15b) 

Equating two relations, we obtain an algebraic equation for w(0). Consequently, u(0) 
is determined by solving the algebraic equation as 

(D - l)n(P) 2 + 2Du(0) + D = =>• u(P) = (-D ± Vl>)/(D - 1), (3.16) 

which implies H 

G 2 (P)-2G(P) + 1-D = P, G(P) = l±v / ^D. (3.19) 



2 The author would like to thank Prof. Taichiro Kugo for correspondence on this issue. 
3 If we do not assume w(0) = 0, we have 



G(0) = 1 + (1 - D)w(0)/2 + VP + (1 - D)w{0)/2} 2 - 1 + D > 0. (3.17) 

by solving 

G 2 (0) - [2 + (1 - L>)w(0)]G(0) + 1 - D = 0. (3.18) 
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Note that u(k 2 ) cannot be smaller than —1, i.e., u(k 2 ) > —1 for G{k 2 ) to be non- 
negative, G(k 2 ) > 0. 

For D = 4, 3u(0) 2 + 8u(0) + 4 = has solutions u(0) = -2/3 and -2. We must 
adopt the solution u(0) = —2/3. Thus we have 

u (0) = -- 0.66666..., G(0) = [1 + w(0)]~ 1 = 3, (3.20) 

irrespective of the number of color N. Thus the ghost dressing function G(k 2 ) is finite 
even in the deep infrared limit k 2 — > 0. This value seems to agree with the bare value 
obtained in numerical simulations on a lattice, see Figure 5 of [12] for SU(2) and 
Figure 4 of [13] for SU(3) (and Figure 2 of [H] for the renormalized values). The 
renormalization point dependence will be discussed later. 

For D = 3, 2u (0) 2 + 6m (0) + 3 = has solutions u(0) = (-3 ± v / 3)/2. We must 
adopt the solution 

u(0) = (-3 + V3)/2 0.63397..., G(0) = [1 + w(O)]" 1 = 1 + V% = 2.73205... 

(3.21) 

This should be compared with numerical simulations on a lattice, e.g., [15], where the 
renormalization must be properly taken into account (separately). 

For D = 2, the derived relations are also valid. However, in two-dimensional field 
theory, a subtle problem exists. As is known in the Coleman theorem for the absence 
of the Nambu-Goldstone particle in two dimensions, a massless particle cannot be 
well defined for D = 2. If the ghost dressing function for D = 2 is finite, this remark 
will be applied. In fact, the numerical simulations indicate that the ghost dressing 
function is not finite and the ghost propagator becomes more singular than the free 
one for D = 2 (HE)]. 

The obtained result differs from the old Gribov prediction in the Landau gauged 
It is possible to reconcile this result with the original Gribov argument as follows. In 
order to see where the difference comes from, we expand the horizon function 

/ 1 \ CE 

yh{x) := 7 / d D ygf ABC ^{x) _ Q (x, y)gf AFE ^ F (y) 



--1 



J d D ygf ABC ^ B {x) (— L-) 5 CE 5^(x - y)gJ AFE ^{y) + 0{{gAf) 

=HNg*?{x) (— + 0((gAf). 

Then the quadratic part in the gauge field of the Gribov-Zwanziger Lagrangian is 
modified as 



(-d 2 ) + 7 i\y 



Jg 2 ; -? + 0((gAy 



The resulting gluon propagator up to 0((gA) 2 ) vanishes in the IR limit: 



-d 2 ) + 7 AV 



1 



-d 2 



-d 2 k 2 



(-d 2 ) 2 + Ng 2 i (k 2 ) 2 + Ng 2 -f 



I (A; i 0). 



4 In the Coulomb gauge, the situation is different and there are at present no contradictions with 
the Gribov prediction in the Coulomb gauge [5T] , 
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In this way the Gribov prediction is reproduced from the Gribov-Zwanziger theory. 
However, we can raise the question: Do the higher order terms modify the IR behavior 
from the Gribov prediction? To answer this question, the full or exact estimation of 
the horizon function is needed as we have done in this paper. 

In our result, the Gribov result was obtained by taking into account only the 
lowest order term as shown in the above. The formal power series expansion, @ 
1+ ^( ) — [1 + ^O)]" 1 = 1 ~~ M (0) + w(0) 2 + • • • yields the horizon condition: 

(h(0)) = (N 2 - 1) {-Du(0) + u(0) 2 - u(0) 3 + ■■•} = (N 2 - l)D. (3.22) 

If we took into account only a linear term in u(0) = 0(g 2 ), then the horizon condition 
would lead to the Kugo-Ojima criterion u(0) = —1 and the divergent ghost dressing 
function G(0) = [1 + m(O)]" 1 = oo. In this way we can reproduce the Gribov original 
result in the Landau gauge in the lowest order. This result totally changes if we 
include the effect of higher order terms. Thus, the ghost propagator behaves like 
free at low momenta. Whereas the gluon propagator is expected to be non-vanishing 
at low momenta where the gluon dressing function vanishes in the IR limit. This is 
possible as shown in [26] . 



4 Proof 

The proof of (13. 2p given in [7J in the case of 7 = is extended to the case of 7 7^ as 
follows. Although the relations we establish in the following have the same forms as 
those given in section 2 of [7J, the Green (or correlation) functions have the implicit 
7 dependence, and hence the average should be understood to be 7 dependent (...) 7 
in the Gribov-Zwanziger theory, which is however omitted for simplicity. 



4.1 The relation (2.5) 

The same relation as (2.5) of [7J in the case of 7 7^ 0: 

((D^) A ^ B ) k = i^5 AB 

is derived in the same way as in the 7 = case. 

In the path-integral quantization, it is derived as follows. 

(0\d»(D^) A (x)tf B (y)\0) 

dti{<5>)e- SGZ d»{D^) A {x)tf B {y) 

5e~ Saz 



<i/i($) 



-S AB 6 D \ 



5(-tf A {x)) 
x-y), 



e- Saz< rf B (y) 



d{ie 



-SrJV B (y) 



(4-1) 



(4.2) 



5 It should be noted that this expansion is not the same as the perturbative expansion in the 
coupling constant. Therefore, this result does not necessarily agree with the result of loop expansions. 
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where we have used a fact that the integration of the derivative is identically zero: H 

Another derivation is given in Appendix B. 
4.2 The 1PI part 

The localized Gribov-Zwanziger theory has the 6 full propagators: 

BC ry/A FBC 



cAB _ FEF 
S^t Si/ ? 

,AB -EF 



ur-ar. ( 4 - 4 ) 



and 5 full vertices: 



<-«f -Of', (4-5) 

where in the tree level only the propagators {^ A {x)^ BC \y)) and (^ A (x)^ BC (y)) have 
the 7 dependence. We denote a set of all fields by $ := ^, £, £, u;, a)}. 

We consider the definition (2.6) and (2.7) for the 1PI part of ((g^ x ^) A< ^ B ) fc 
and l ^€ A {gsrf v x c £) B ) k . For the 3-point function of gluon, ghost and antighost, there 
is no disconnected part. See Fig. [T] for a diagrammatic representation. In defining 
the 1PI part, the possible intermediate fields $!,$ 3 ,$ 5 are uniquely determined to 
be $3 = ff, $ 5 = ^ and $i = stf by taking into account the propagators and vertices 
enumerated in the above. Therefore, the 1PI part is immediately defined by 

((g^(x) x tf{x)) A( i B {y)) ={{g^(x) x tf(x)) A tf c (z)) lpl (tf c (z)tf B (y)). (4.6) 

In the similar way, we can define 

^ A {x){g^{y) x tf(y)) B ) ={V A (x)V c (z))(<<f c (z)(g^(y) x <f(y)) B ) m . (4.7) 

Then the Fourier transform reads 

((£< x ^) A ^ B ) k =((g^ x tf) A tf c )l P1 (tf c tf B ) k , (4.8) 

(K A {gK x <i) B ) k ={<£ A <£ c ) k {<<? c (gK x ^) B )\ P \ (4.9) 



6 This relation is called the Schwinger-Dyson equation, which is a consequence of J a d<f>^ = 
f(b) — f(a) = where we have used f(b) = f(a) at the boundaries. This relation follows also 
from the shift invariance of the measure, since J d<f)f(4>) = / d(<f> + a)f(<fi + a) = J d<frf((j) + a) = 
J #/(</>) + a J d(j)^ + 0(a 2 ) holds for arbitrary a. 
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c c 

c w c 



1PI 



Figure 1: Diagrammatic representation of (a) ((g^ x c &) Ac & B )k and 
^) A( i B )f\ (b) (<^(<7< x and {^ A (gK x < #) B )^ PI . 

4.3 The identity (2.10) or (2.11) 

An identity (2.10) or (2.11) of [7J: 

= {{d»D^) A {gsrf v x V) B ) k (4.10) 

is necessary to define the Kugo-Ojima function according to (13.51) . This is derived in 
the case of 7 ^ in the same way as in the 7 = case. The transversality follows 
from the identity: 

= / M®)^^ le- SGZ (9< x <£) B {y)\ 

d^)e^[gf ABC K C {y)5 D {x -y)- —j^** * <£) B {y)\ 
= J d^)e^[gf ABC K°(y)S D (x -y)- d»(D^) A (x)(gK x tf) B (y)}, (4.11) 
which leads to 

{^{D^) A {x){gK x V) B (y)\0) = gf ABC (0\^ u c (y)\0)5 D (x - y) = 0, (4.12) 
since (0| I e^p(?/)|0) = from the Lorentz invariance of the vacuum. 

4.4 The relations (2.13)-(2.15) and the 1PI part 

We proceed to check the relations (2.13)-(2.15) of [7J. 

In order to examine the constraint coming from the horizon condition in the 
Gribov-Zwanziger theory, the most non-trivial issue in discriminating between 7 = 
and 7^0 cases is how to connect the function ({gs^u x ^) A (g£/^ x ^) A )fc in the 
average of the horizon function (h(0)) to the 1PI parto 

(h(0)) = V D l j d D x(h(x)) = - lim((^ x <tf) A (g^ x <tf) A ) k . (4.13) 

J k— »0 



7 Here we have used the identity which holds for any functional f(A) of A: 
(f(A)V A (x)V B (y)) = -(f(A)(K-y B (x,y)). 
This is consistent with (|4.2p . 
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ml P I s i ngl 



Figure 2: Diagrammatic representation of ((g^ x ( &) A (g£/ u x c io) B )k, ((g^ x 
^) A (#< x tf) B ) c k onn , ((g^ x x <i) B )f l and ((^ x <&) A {gst v x ^) B )^ 1PI . 

In what follows, we check that no 7-dependent extra terms do appear in this 
process even in the Gribov-Zwanziger theory. As a result, the relations derived in 
the previous paper [7j still hold with no change even in the case of 7 7^ 0. In other 
words, the relations and the results derived from them in the previous paper [7] are 
valid irrespective of the value of 7. 

The 4-point function of gluon, gluon, ghost and antighost is separated into the 
connected part and the disconnected part, see Fig. [21 

((^On) x <£{x 2 )) A {gK{yi) x <% 2 )) B > 
:=gf AFC gf BGE (^ F (x 1 )K G (yi))(^ c (^W E (y2)) 

+ ((jgKM x ^2)) A (gK(yi) x V(y2)) B ) conn . (4.14) 

After putting x\ = X2 '■— x and y\ = y<i := y and taking the Fourier transform, we 
obtain the necessary relations. The connected part is defined by El 

(GK x <g) A {gtf v x <£) B ) k 

■= J ^9f AFC 9f BGE KK G ) k+I ^ E ) P +{{9^ x V) A {gK x <i) B rr- 

(4.15) 

8 Eq.(2.15) in [7] involves misprints in the indices. Here we give correct expressions. 
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The 1PI part is defined from the connected part as shown in Fig. [5] by using a 
diagrammatic representation. Here we use the same type of argument to define the 
1PI part. Consequently, the 1PI part is defined from the connected part 

((gj*p x V)\gK x <i) B Yr = ((9< x tf)\gK x <£) B )™ + A AB (k), (4.16) 

where A AB involves the single ghost propagator: 

Ku( k ) ■= (GK x ^) A ^ C )f\^ C ^ D )k^ D {9< x <g) B Y*\ (4.17) 
Therefore, 1PI and mlPI is related as 
-M(9< x V) A {gK x V) B )f pl 

■= J ^9f AFC 9f BGE (^ F K G h +P (^ E ) P +((9K x x V) B )f 

(4.18) 

Then, we can write the relation 

((^ x syW, x ^) B ), = Ajf (A) + A AB (k). (4.19) 
Thus the average of the horizon function is written as 

(MO)) = - Hm((^ x *) A ^ x V) A ) k = -X AA (0) - A^(0). (4.20) 

4.5 The relations (2.16)-(2.25) 

The relations (2.16)-(2.25) of [7J remain unchanged. 

In the manifestly covariant gauge of the Lorenz type, an idenity holds ((2.17) 
of 0) 

ik^{k) = (V A (g< x <i) B )f\ (4.21) 
In the Landau gauge, a similar identity also hold ((4.1) of [7]) 

A£?(A0(-ik) = (( 9 K x V)***)? 1 - ( 4 -22) 

This identity (14.221) is derived from (14.21j) using the Faddeev- Popov conjugation in- 
variance in the Landau gauge. See Appendix A. 

In addition, an identity holds in the manifestly covariant gauge, ((2.3) or (2.8) 
of 0) 

((g^ x <£) A <i B )™ = -ik, (-8 AB + -±(V*f")j;^ , (4.23) 

yielding 

ik^g^ x tf) A tf B )l P1 = -5 AB k 2 - OrtP)* 1 , (4.24a) 
which is the same as the Schwinger-Dyson equation for the ghost propagator ((2.9) 

of my- 

1 un 

(V A V B ) k = -5 AB - - i^{(g^ x V) A V B ) k . (4.24b) 
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Using the general form for the uncontracted 2-point function X AB (k) ((2.25) of [7J 



1 1 V 



5,Mk 2 ) + ^w(k 2 ) 



irAB 



(4.25) 



Z,' 2 

we obtain a relationship between the Kugo-Ojima function and the ghost propagator: 
G(k 2 )5 AB := -k 2 i^ A( i B ) k = [1 + u(k 2 ) + ^(A; 2 )]" 1 ^. (4.26) 
Incidentally, the relation (14.251) is obtained from (13.51) and ((2.20) of [7]) 



((D^) A (g*f v x tf) B ) k = ( 5 W - ^ ) (GK x ^)\g^ v x tf) B )f PI . (4.27) 



Here w(k 2 ) is an unknown function. By combining (14.251) with (I4.2ip . (I4.22p and 



04.261) . we obtain 



=(iV 2 - 1)[/Mfc 2 ) + w(fc 2 )], (4.28) 



= - (iV 2 - l)G(A; 2 )[u(A; 2 ) + w(A; 2 )] 2 

v 7 1 + u(fc 2 ) + «;(&; 2 ) v y 

The existence of the last term A AA (0) is crucial to obtain a finite ghost dressing 
function, see 06 .81). 



5 Renormalization 

In order to compare our result with the numerical simulations on a lattice and the 
Schwinger-Dyson equations, we need to consider the renormalization, especially, the 
dependence of the renormalization point under a given renormalization condition. 
We denote the gluon propagator in the Landau gauge by 

rtAB(i\ x ab I x k^k u \ F(k 2 ) 

D^{k) = S l<W--p-J-p-> I 5 - 1 ) 



and the ghost propagator by 

a AB (h\ = s AB - 

k 2 ' 



G AB {k) = -8 AB ^p-. (5.2) 



The gluon-ghost-antighost vertex is denoted by 

T ABC (p,k) = f ABC T u (p,k). (5.3) 
In the bare case, d(k 2 ) = 1/k 2 , G(k 2 ) = l/k 2 and T u {p,k) = —k u . 
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By using the identity (13. SJ) , 



BxmlPI 



S^ u u(k 2 



^ W (k 2 ) 



S 



AB 



(5.4) 



it is easy to show that the functions u(k 2 ) and w(k 2 ) obey the equations: 

1 



u(k 2 ) 
w(k 2 ) 



(D-1)(N 2 
1 

(D-1){N 2 



A^(fc) - ^A^(fc) 



ht JL h u 

ft fi, AA 



k 2 
(k) 



Km 



(5.5a) 
(5.5b) 



The momentum dependence of the functions u(k 2 ) and w(k 2 ) is determined from the 
knowledge of the function: \^{k) := ((#^x^) A (#<x^) B )™ lpI , see Fig.[S This is 
performed e.g., by substituting the solutions of the Schwinger- Dyson equation for the 
propagators and the vertex functions into the right-hand sides of (I5.5al) and f!5.5bl) . 
as done in [23] . 

We consider the multiplicative renormalization in which the renormalization con- 
stants are introduced in the standard way: 



F R (k 2 ^ 2 ) =^V,A 2 )F( K 2 ,A 2 ) 
G R (k 2 ,fi 2 ) =Zc\n 2 ,A 2 )G(k 2 ,A 2 ) 
r» R (k,P,H 2 ) =Z 1 ( yU 2 ,A 2 )r(A;,p,A 2 ) 
g R (ji 2 ) =Z~ 1 ( f i 2 ,A 2 )g( f i 2 ), 



(5.6a) 
(5.6b) 
(5.6c) 
(5.6d) 

where /i is the renormalization point, A is an ultraviolet cutoff, and Z g = Z ^^ 2 Zq X Z\. 

We carry out the renormalization in such a way as to preserve the identity (13.101) . 
In order to preserve the identity (I3.10p after renormalization, we impose that 



G R \k 2 ^ 2 ) =Z c (^ 2 ,A 2 )G~ 1 (k 2 ,A 2 ), 
1 + u R (k 2 , fi 2 ) =Z c (fi 2 , A 2 ) [1 + u(k 2 , A 2 )], 
w R {k 2 ^ 2 ) =Z c ^ 2 ,A 2 )w(k 2 ,A 2 ), 



(5.7a) 
(5.7b) 
(5.7c) 



Consequently, the Kugo-Ojima function u R (k 2 ) is not a renormalization-group 
invariant quantity. Hence, it is not a renormalization point /i-independent quantity 
and the Kugo-Ojima parameter u R (0) depends on the renormalization point /i. In 
the case of an IR finite ghost dressing function, therefore, u R (0) acquires a non- 
trivial dependence on the renormalization scale, while in the case of an IR divergent 
ghost dressing function the possible /i-dependence is not important, since u R (0) = — 1 
irrespective of the value of /i chosen. 

Moreover, we introduce the renormalization constant Z\ which relates the bare 
and renormalized functions, A^f and X^ R , as 



^JAk,^ 2 ) = Z x ^ 2 ,A 2 )X A J(k,A 2 ) 



(5. 
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We find that by virtue of the identities which are the same as (2.17) and (4.1) 
of EZJ: B 

W\fi(k) =ik»((g^ x V) A (gK x tf) B )f PI = (V A (gK x <f) B )™ \ (5.9a) 
-iAif (*)*" = - i((g^ x Vfitj*, x <*f) B )f PI r = ((^ x V)***)™, (5.9b) 



in the Landau gauge X AB (k) must be renormalized by the same renormalization con- 
stant as that for (tf A (g^ x <£) B )f l (and ((g^ x <T) A< # B )£ PI ) which appears in the 
Schwinger-Dyson equation for the ghost dressing function or the ghost propagator: 

G AB (k)~ l := {^ A( i B )l l = -k 2 5 AB - ik"({gsfp x ^) A< tf B )f l . (5.10) 

Therefore, we have 

Z A (/i 2 ,A 2 ) = Z (Ai a , A 2 )£rV, A 2 ). (5.11) 
The renormalized quantities must obey 



l + u R {k\^)=Z c {^\^ ' Zl(/i2 ' A ' 



(D - l)(iV 2 - 1) 



(5.12a) 



r ,2 2 ^ ^i(/i 2 ,A 2 ) 



(£> - l)(iV 2 - 1) 



(5.12b) 



<#(*V) =Z c (fJ 2 ,A 2 ) + ^^ ^lA^(fc )A i a ). (5.12c) 

If Zc is determined, we can cast the SD equation for Gn(k 2 ) := k 2 Gn(k 2 ) into a 
manifestly renormalizable form. The same Zq makes un(k 2 ) finite. 

Once the renormalization condition for Gr is chosen to be Gr(/i 2 ) = 1, the values 
of uji(fi 2 ) is completely determined for its own equation (15.1 2ap . It should be noted 
that the horizon function (Air(0)) is not a renormalization-group invariant quantity, 
and hence depends on the renormalization point just as the Kugo-Ojima parameter 
ur(0) depends on the renormalization point \i. 

We now study the //-dependence of (the average of) the horizon function. 

(h(0)) = -\ AA (0)-A AA (0). (5.13) 



If the renormalization of A AB (k) := ((gK x ^) Ac € c )f l ^€ Cc € D ) k ^ D {g^ v x ^) B ) 



ipi 



is done according to A AB R {k) = Z&A AB {k) as 

KU k ) = (ZiZ^Z^Z^Z^Z^A^k) = Z^Z c A AB (k) := Z A A AB (k), 

(5.14) 

then the horizon function is renormalized as 

(MO)) = -A^(0) - A AA (0) = -Zc\Z x \ AA R ($) + Z 2 A AA R (0)]. (5.15) 



9 The relation used in |23j k^H^ ll ,{p 1 k) = —iT^(jp,k) must be modified in the Gribov-Zwanziger 
theory. 
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The Taylor's non-renormalization theorem [32] says that in the Landau gauge, the 
renormalization constant Z\ of the ghost-antighost-gluon vertex is exactly one, Z\ = 
1, see also section 2.2 of [33] . 

The /z-dependence of the horizon condition would be obtained from the renormal- 
ized horizon condition: 

(MO)) = -(N 2 - 1) {Du R (0) + w R (0) - G R (0)K(0) + w R {0)} 2 } . (5.16) 

This should be compared with the renormalized expression which comes from another 
form of the bare horizon condition given in a previous paper [7], 

(MO)) = -(N 2 - 1) {{D - 1)M0) + G R (0)[u R {0) + w R (0)}} . (5.17) 

The issue of dependence of the Kugo-Ojima function and the horizon function 
on the renormalization point \i has been studied in detail by Aguilar, Binosi and 
Papavassiliou [23] using a mixed approach, in which the lattice data for the gluon 
propagator were used as input for their SD equations, see Fig. 9. According to their 
results, the renormalized values which are the same as the bare value u(0) = —0.665 ~ 
— 2/3 (and (h(0)) = 31.92 ~ (3 2 — 1)4 = 32) are obtained when the renormalization 
point is chosen at ji ~ VWmGeV ~ 3.3GeV. It is possible to translate this result 
to that at the other renormalization point by using their result. 

The equation for w R (k 2 ) remains unchanged, i.e., one simply replaces in the equa- 
tion for w(k 2 ) the unrenormalized quantities with the renormalized ones, up to Z\. 
They have shown that w R (0) = if both the gluon propagator and the ghost dressing 
function are IR finite under the approximation that the form factors are replaced with 
their tree level values, i.e., bare vertex approximation. The assumption appears to 
be a good approximation according to lattice studies [23130] • Note that the pertur- 
bation theory leads to w(0) ^ 0, e.g., w(k 2 ) = g 2 N/ (32tt 2 ) at one loop, because the 
free gluon propagator is not IR finite. 

The renormalization will be discussed in more detail in a subsequent paper [38J. 

6 Remarks 

6.1 Remark 1: The total derivative term in a choice of the 
horizon term 

We point out that the result crucially depends on the explicit form of the non-local 
horizon term adopted. 

If the total derivative was neglected in the Gribov-Zwanziger theory, the horizon 
term could be rewritten as 

/ d D xh(x) := J d D x J d D ygf ABC ^{x){K-r E {x,y)9f AFE ^{y) 

?=Jd D xj d D yD^} AC (x)(K~ 1 f E (^y)9f AFE K F (y) 
^=Jd D xJ d D yD^] AC (x)(K- 1 ) CE (x,y)D tM [^] AE (y), (6.1) 

10 The author would like to than Dr. Daniele Binosi for the correspondence on this point. 
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which yielded the average of the horizon function: 

(h(0)) = - hm((D^) A (D^) A ) k , (6.2) 

instead of (h(0)) = — lim^oU^- 2 ^ x ^) A (fi'^ x ^) A )k- Thus, the horizon condition 
(h(0)) = (N 2 — 1)D for this horizon function defined from this form, i.e., 

(MO)) = -Km((D^) A (D^) A ) k = -(iV 2 - 1) {(D - 1)«(0) - 1} = (iV 2 - 1)D, 

(6.3) 

led to the Kugo-Ojima criterion and the divergent ghost dressing function: 

u(0) = -1, (6.4) 

and the divergent ghost dressing function with an imput w(0) =0: 

G(0) = [1 + u(0) + w(O)]" 1 = ^(O)' 1 = oo. (6.5) 

Thus, if one starts from the horizon term in the last form of ( 16.11) as done in [3lj 
by neglecting the total derivative term in the non-local horizon function, then one 
is led to the opposite conclusion to ours. This fact that the horizon condition using 
this definition (16.31) is equivalent to the Kugo-Ojima criterion (16. 4p has already been 
pointed out and it was checked to what extent the horizon condition holds in the 
numerical simulation in [33] • More details on the foundation for the horizon function 
are given in [3"o] . 

The resulting discrepancy comes from the neglection of the total derivative term 
in the non-local theory. In the case of the local theory, the total derivative term can 
be neglected, since the volume integral is replaced by the surface integral over the 
surface S far away from the origin by the Gauss theorem: 

/ d D xdMx)] = [ d^SpWlfpiz)] 

JV JS=dV 

= Jdn D \x\ D - 1 nMx)=0, (6.6) 

for a rapidly decreasing and vanishing function at infinity: n^f^x) ~ \x\~ a (a > 
D — 1), where flo is the solid angle. 

In the case of the non-local theory, however, the total derivative term cannot be 
neglected, even if the integrand f(x — y) is a rapidly decreasing function in \x — y\: 

f d D y f d D xd;[f(x-y) gix (y)\ = f d D y f d D ~ l S,{x)[f{x - y)g,{y)} 

JV JV JV J S=dV 

= [ _ d^S^x) [ d D y[f(x - y)g,{y)] ± 0, 

(6.7) 

because a contribution from y such that [re — j/j ~ always remain due to the existence 
of an additional integral J v dPy over the whole volume V. Even after applying the 

11 The author would like to thank Drs. David Dudal and Nele Vandcrsickel for kind discussions 
on this issue. 
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Gauss theorem, we have still another volume integral and the integrand is not limited 
to the region far from the origin. Thus, if one neglects the total derivative term in 
the non-local theory, one is led to a theory with different physical content. 

It is a good place to note that, in order to obtain the modified BRST symmetry 
in [2T1I22] . integration by parts, i.e., neglection of the total derivative has been used 
as usual. However, this is safe, because it is done in the localized Gribov-Zwanziger 
theory. 

6.2 Remark 2: A contribution from a diagram connected by 
a single ghost line 

It should be remarked that our result differs from the Zwanziger result [37]. The es- 
sential difference between his result and ours comes from the estimation of A sing (k) := 
Avjf (k). In fact, if the contribution from this term was neglected, the average of the 
horizon function became equal to (instead of (I3.15al) ) 

(MO)) = -A#(0) = -(N 2 - l)Du(0), (6.8) 

and the horizon condition (h(0)) = (N 2 — 1)D would become indeed equivalent to 
the Kugo-Ojima criterion 

u(0) = -1, (6.9) 

as claimed in [37]. However, A sing (k) does not vanish and remains non-zero even after 
taking the k = limit, as we have examined in the previous paper [7J, which is also 
the case for 7 ^ 0. 
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A Symmetry of exchanging ghost and antighost 

In the Landau gauge, d^ A (x) = 0, the equations of motion for the ghost field and 
the antighost field have the same form 

d^D^\€ A = 0, d^D^\€ A = 0, (A.l) 
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This suggests that there is a symmetry under the exchange of ghost and antighost 
fields. In fact, the GZ action is invariant under the FP conjugate transformation: 

where other fields are unchanged and SB is defined by 

S8 A + m A = g(itf x ^) A . (A.3) 
This is easily seen by rewriting the GF+FP term as 

Sgf+fp := / d D x [@ ■ 8^ + • d^Dtf} := J al D x [-8 ■ - V ■ d»Dj<£} . 

(A.4) 



B Another derivation of eq. (14.21) 



We work in the Minkowski formulation. In the canonical quantization, the relation 
g2D, i.e., 

d%{0\T(D^) A (x)V B (y)\0) = 6 AB 6 D (x-y), (B.l) 

is derived from the equation of motion for the ghost field d^D^ix) = and the equal- 
time canonical anti-commutation relation{^ A (x), i(Do<o) B {%))} Xo=yo = i8 AB 8^ D ~ X > {x— 
y). In fact, the equation of motion 

and the equal-time canonical anti-commutation relation 

i5 AB 8 D -\x -y) = {Tl A -(x) : V B (y)} Xo=yo , (B.3) 

implies 

dZ(0\T(D^) A (x)V B (y)\0) 
=5(x - y )(0\{(D V) A (x)tf B (y)}\0} + (0\Td»(D^) A (x)V B (y)\0) 

=5 AB 5 D (x-y), (B.4) 

where we have used 

r) 9^ 

^ = mm = i(D ^ )A{x) ' (R5) 

and the time-ordered product (T-product) for the anti-commuting field defined by 
TA(x)B(y) := 6(x - y )A(x)B(y) - % - x )B(y)A(x). (B.6) 
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